mt<-read.table("tmc_mass.txt",head=FALSE)
mdt<-density(mt[[1]])
mdt[[2]]<-mdt[[2]]*18
plot(mdt, xlim=c(0,10),xlab="Mass/M\u2299", ylab="dN/dM", main="CMF of TMC")
abline(2.8,-0.3,lty=2)
abline(v=1,lty=3)

mt<-read.table("cmc_mass.txt",head=FALSE)
mdc<-density(mt[[1]])
mdc[[2]]<-mdc[[2]]*14
plot(mdc, xlim=c(0,100),xlab="Mass/M\u2299", ylab="dN/dM", main="CMF of CMC")
abline(0.2784,-0.00234,lty=2)
abline(v=40,lty=3)





























































































































































































































































































































































































































































































































